Astronomy Reports, 2010, Vol. 54, No. 4, pp. 338—354. 



The Appearance of a Radio-Pulsar Magnetosphere 
from a Vacuum with a Strong Magnetic Field. 
Motion of Charged Particles 

Ya. N. Istomin^ and D. N. Sob'yanin^ 

^Lebedev Physical Institute, Russian Academy of Sciences, 
53 Leninskii pr., Moscow, Russia 

"^Moscow Institute of Physics and Technology (State University), 
Dolgoprudnyi, Moscow oblast, Russia 

Received June 30, 2009; in final form, October 19, 2009 

Abstract 

Tlie motion of electrons and positrons in the vacuum magnetosphere of a neutron 
star with a surface magnetic field of i? ~ 10^^ G is considered. Particles created in the 
magnetosphere or falling into it from outside are virtually instantaneously accelerated 
to Lorentz factors 7 ~ 10^. After crossing the force-free surface, where the projection 
of the electric field onto the magnetic field vanishes, a particle begins to undergo ultra- 
relativistic oscillations. The particle experiences a regular drift along the force-free surface 
simultaneous with this oscillatory motion. 
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1. INTRODUCTION 



The stationary structure of the magnetosphere of a radio pulsar has been studied in 
considerable detail. Here, this refers not to the vacuum magnetosphere, but instead to a 
magnetosphere filled with dense electron-positron plasma. This is associated with the fact 
that the radio emission generated in the magnetosphere by the flux of charged particles 
requires the action of some mechanism producing a dense plasma. In one such mechanism, 
proposed by Sturrock [l!] and substantially developed by Ruderman and Sutherland 
the efficient creation of electron-positron pairs by gamma-rays with energies exceeding 
twice the electron rest-mass energy is possible in the strong magnetic field at the surface 
of the neutron star — radio pulsar, B ~ 10^^ G [s^, !4|. In turn, energetic gamma-rays 
are emitted by electrons and positrons during their motion in the magnetosphere along 
magnetic-field lines having appreciable curvature. Such photons are called curvature 
photons. The emission and absorption of photons in a magnetic field accompanied by the 
creation of electron-positron pairs provides a theoretical basis for the steady-state creation 
of plasma in the magnetosphere of a magnetized, rotating neutron star. The density 
of the forming plasma n appreciably exceeds the so-called Goldreich- Julian density |5|, 
ugj = |f2 • B|/27rce, which provides stationary rotation of the magnetosphere right out 
to the light cylinder, Rl = c/Q. Here, Q is the angular velocity of the star's rotation, 
c is the speed of light, and e is the positron charge. The multiplicity parameter for the 
creation of the electron-positron plasma A = n/ncj is large: A ~ 10^—10^. 

However, it is difficult to understand from observations of stationary radio pulsars 
what radio-emission mechanisms are operating in their magnetospheres, and where and 
how the plasma is produced. Tracing the dynamics of the development of the emission at 
different frequencies could be very important for our understanding of the physical pro- 
cesses occurring in the magnetospheres of radio pulsars [g]. Moreover, many observations 
of non-stationary radio pulsars have recently appeared. These are first and foremost so- 
called switching radio pulsars, from which radio emission is observed only during certain 
time intervals appreciably exceeding the rotational period of the star. For example, PSR 
B1931+24 is "on" for 5-10 days, then "off" for 20-25 days fzi; PSR J1832+0029 is "on" 
for about 300 days and then "off" for about 700 days (see, e.g., the review [sl). Measure- 
ments indicate that the deceleration of the rotation, i.e., the loss of energy, is appreciably 
different during the "on" and "off" periods. Since the power of the radio emission is a 
negligible fraction of the total rotational energy lost by the neutron star, it is natural to 
suppose that quiet periods are associated with the cessation of the generation of plasma in 
the magnetosphere. We can then distinguish loss mechanisms associated with the radia- 
tion of magnetodipole waves in vacuum ( "off" periods) and with the emission of the pulsar 
wind in the form of a flux of electron-positron plasma flowing from the magnetosphere 
("on" periods). 

In addition to switching pulsars, a group of so-called nulling pulsars has long been 
known, which likewise display no radio emission during certain intervals, but with these 
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being not as regular as those for switching pulsars. Differences in the rotational decel- 
eration have not yet been measured for nulling pulsars. For example, PSR B1944+17 
does not have detectable radio emission 55% of the time 0. Of the 23 pulsars studied in 



10| . 7 display nulling fractions exceeding 40%, while this fraction reaches 95% for PSR 



J1502-5653 and PSR J171 7-4054. 

Another group of non-stationary radio sources has recently been observed: rotating 
radio transients (RRATs), which are sporadically flaring radio sources. The phases are 
preserved during these flares, and the corresponding measured periods are characteristic 
of ordinary radio pulsars [11]. There is no doubt that these are also rotating neutron stars. 
However, the nature of their activity is quite unclear, as is testified to by the presence 
of numerous and varied models for these objects, such as models invoking precession 

"dead" pulsars 



, re-activation o 



12l. reversal of the direction of the radio beams 113 

13|. interaction of the magnetosphere with a disk |l4|, drift waves [l5|, and even such 
exotic objects as the remnants of quark novae [1^. Naturally, such models are open for 
discussion. For example, when considering re-activation of the radio emission of pulsars 
located to the right of the "death line" in the P—P diagram, we must know the exact 
position of this line. Timing data for the RRAT J18 19-1458 suggest that the magnetic 
field at the stellar surface is approximately 5 x 10^^ G, so that it exceeds the Schwinger 
field jll, 17|. In such a strong magnetic field, it is necessary to consider the splitting of 
photons taking into account their polarization and different conditions for the creation 
of pairs than in the case of a weak magnetic field. As a consequence, the death line for 



pulsars with strong magnetic fields has a slope of 11/3 18|, ll9j rather than the usual slope 
of 11/4, which occurs only in the case of weak magnetic fields. 

In our view, all these sources exhibit a non-stationary generation of plasma in 
neutron-star magnetospheres. It is therefore important to understand how the magneto- 
sphere of a rotating, magnetized neutron star is filled with plasma, which is the reason 
for the operation of radio pulsars. An important aspect of this is the state of the magne- 
tosphere of the rotating neutron star before it becomes a radio pulsar, or after the source 
of plasma has switched off in its magnetosphere. If this source does not operate, the 
dense plasma falls onto the surface of the neutron star in a closed magnetosphere over the 
characteristic time t ~ Rl/c = = P/2'k, which is less than the rotational period. In 
an open magnetosphere, the plasma flows out over the same time scale. This raises the 
question of what then remains in the magnetosphere. 

One important discovery here was yielded by observations of the dynamics of the 
braking of the two switching pulsars PSR B1931+24 and PSR J1832+0029 0, fl, which 
showed that the deceleration of the star's rotation continued after the disappearance of the 
radio emission. The rate of this deceleration was below the initial value by approximately 
a factor of 1.5, but not by an order of magnitude. It is not possible to explain the observed 
braking if plasma with a density of the order of ncj remained in the magnetosphere after 
the pulsar switched off, so that the screening of the longitudinal electric field continued. In 



3 



this case, the magnetodipole radiation would be fully screened. If there were energy losses 
associated with the outflow of plasma with the density ugj from an open magnetosphere, 
the particle energies would have to be of the order of 10^^ eV, which would require the 
presence of an acceleration region in which the longitudinal electric field was not equal to 
zero. An energetic production of plasma with multiplicity A ^ 1 would then begin in this 
region, in contradiction with the fact that the pulsar was switched off. Moreover, it may be 
that there is no plasma in an open magnetosphere, but that it occupies some closed region. 
However, as was shown in [gI, magnetodipole radiation will be weakened by a factor of 
{R/RlY/'^ ~ 10-6-10-^ < 1 in this case, compared to the case of a pure vacuum, which 
is not observed. The only reasonable conclusion is that there is no plasma in the pulsar's 
magnetosphere when it is not in its operational state, and that the observed rotational 
deceleration is brought about via magnetodipole radiation. In this case, the energy losses 
have the same order of magnitude as for the operational pulsar, when the generation of 
plasma and electric current in the magnetosphere lead to braking of the rotation. We do 
not discuss here what this would imply for various plasma-production models developed 
over many years; the absence of plasma in the switched-off state (n <^ ^gj) represents 
direct evidence from available observational data, and not a new proposed model. 

We are concerned here with an initial investigation into the "ignition" of the mag- 
netospheres of neutron stars — the dynamics of the filling of the vacuum magnetosphere 
with electrons and positrons created in the magnetosphere. The second section considers 
the electromagnetic fields of the inner vacuum magnetosphere and describes the force-free 
surface where the electric field component along the magnetic field vanishes. The follow- 
ing sections are dedicated to the dynamics of the motion of charged particles near the 
force-free surface. 



2. STRUCTURE OF THE VACUUM MAGNETOSPHERE 

It will be convenient for us to consider the electromagnetic field around a neutron star in 
the spherical coordinates (r, 6', yj). We take the polar axis to coincide with the rotational 
axis of the star, which determines the direction of its angular-momentum vector f2. Here, 
r is the distance from the center of the star to a given point, 9 the polar angle measured 
from the rotational axis, and the azimuthal angle. The electric and magnetic fields will 
be written in the form 

E = Ere.r + EQe.0 + E^e^, 

B = BrCr + Bgeg + B^B^, 

where the unit vectors 6^,69, are mutually orthogonal and form a right-handed set. 
The electromagnetic field outside the neutron star was found by Deutsch [20(] and has the 
form 
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magnetic field 



2m 

m 

-De — ^ 

m 



cos cos 6*^ + sin 6 sin 6^^ cos(v9 — iprn) 
sin cos — cos sin 6*™ cos(v9 — iprn) 



(1) 



5^ = — sin6'm sin(v9 - Lpm) 



and electric field 



Er 



Eo = -k 



cos 29 + ^ j cos 6'm + ^ sin 29 sin 6'm cos(v9 - iprn) 

2 

sin 26' cos 9m + — cos 29^ sin 6'^ cos(v3 — 9?m 



(2) 



= ~ ^ ) cos 6' sin 6*^ sin (v? - ipm)- 



Here, 



, (fm) are the polar and azimuthal angles of the magnetic axis, which is deter- 
mined by the direction of the magnetic dipole moment m, k = Q/c is the wave number 
corresponding to the angular frequency of rotation of the neutron star Q, and R is the 
radius of the star. The azimuthal angle ipm = is chosen to that it is equal to zero 
at time t = 0. Although we are primarily interested here in the case of an uncharged, 
rotating, magnetized neutron star possessing a high conductivity, the generalized solution 
of Deutsch for a charged sphere can be found in [2l|, |2^ . A useful discussion of this 



question is also presented in the review 23|. The solution of Deutsch refers to a dipo- 
lar magnetic-field distribution at the stellar surface. Its generalization to the case of an 
arbitrary, axially symmetrical magnetic-field distribution is presented in 24 1. 

Formulas ([T]) and ([2j) are valid until we consider the electromagnetic field at distances 
that are appreciably less than the radius of the light cylinder, Rl = c/Q. More precisely. 



(3) 



where r± is the distance from the rotational axis of the neutron star to the point con- 
sidered. As r±_ is increased in ([1]) and (|2]), terms of the expansion with even powers 
of (fir^/c) begin to appear. Further, we will be interested in the generation of elec- 
tron-positron plasma in the neutron-star magnetosphere. As we know, the single-photon 
creation of an electron-positron pair is efficient only in the presence of a fairly strong 
magnetic field. The magnetic field, which has a dipolar structure, rapidly falls off with 
distance r from the center of the star (oc r~^). The creation of electron-positron pairs 
becomes inefficient at magnetic-field strengths of ~ 10^ G. For pulsars with character- 
istic surface magnetic fields 10^^ G, the characteristic distance at which the creation of 
pairs becomes possible is of the order of (10—20)/?. Therefore, for typical stellar radii 
i? ^ 10 km and rotational periods P ~ 0.1— 1 s, the parameter is of order 10^^—10^^. 
Under these conditions, ([1]) and (|2])yield fairly precise results. An impression of the struc- 
ture of the vacuum magnetosphere at distances comparable with Rl can be obtained, for 
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example, from 25 1 



To further analyze the dynamics of the particle motions in the vacuum magneto- 
sphere of the neutron star, we will use the equation for the so-called force-free surface — 
the surface satisfying the equation E ■ B = 0. We are interested in this surface for the 
following reason. We must determine how the charged electrons and positrons that are 
created in the magnetosphere will move in the magnetosphere, and whether there will 
exist regions where they accumulate. We can immediately qualitatively assert (we will 
consider this question quantitatively below) that, after their creation, particles will virtu- 
ally instantaneously be accelerated to relativistic speeds if the longitudinal electric field 
£"11 = E ■ B/i? differs from zero. The only region where particles could collect is the 
force-free surface, at each point of which the longitudinal electric field is zero. 

Using ([T]) and ([2]) for the Deutsch electromagnetic field, it is straightforward to obtain 
an expression for the scalar product of the electric and magnetic vectors (see also [2^): 



2 

m 



E - B = -kr— 



^•6 



( 1 T ) cos 9" sin Om + ^4 cos^ 9' cos 9 



(4) 



Here, we have introduced the angles 9' and 9" as follows: 

cos 9' = cos 9 cos 9m + sin 9 sin 9m cos{(p — (pm), 
cos 9" = — cos 9 sin 9m + sin 9 cos 9m cos((y9 — ipm)- 

Let us define these angles less formally. Let = r/r and = m/m be unit vectors 
directed along the radius vector r and along the magnetic axis m. Further, we define the 
vector e„ = dem/d9m, lying in the plane of the vectors fl and m, directed orthogonal to 
m and indicating the direction of growth in 9m- Then, the angle 9' is the angle between 
the radius vector and the magnetic axis, and 9" is the angle between the radius vector 
and the e„ axis; i.e., cos 6'' = ■ and cos9" = ■ e„,. 

Using (jl]), we can easily obtain an equation for the force-free surface where E-B = 0: 



COS 9 cos^ 9' 
~ V" ~"sin^^ cos^" 

The structure of the force-free surface ([5j) is illustrated by the figure, which shows 
cross sections of the force-free surface by the plane ip — ipm = {0,vr}, passing through 
the rotational axis and magnetic axis, and by the plane ip — ipm = {7r/2, 37r/2}, passing 
through the rotational axis orthogonal to the previous plane, for various angles between 
the magnetic and rotational axes of the neutron star. An understanding of the form of 
the force-free surface can be obtained from j26|. 

In the case of a co-axial rotator, the force-free surface is simply the equatorial plane, 
9 = 9' = 7[/2. For a rotator with an arbitrary inclination, the force-free surface can 
be divided into two regions. One region has the form of two arched parts adjacent to 
the surface of the neutron star at the points of the equator and magnetic equator. The 
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magnetic equator is the circle that is the cross section of the 9' = 7r/2 plane passing 
through the center of the neutron star orthogonal to the magnetic axis, m, and the stellar 
surface, r = R. The ordinary equator is the cross section of the 6 = 7i/2 plane passing 
through the center of the star orthogonal to the rotational axis, fl, and the stellar surface. 
The force-free surface is adjacent to the surface of the neutron star at all points of the 
magnetic equator, which is not true at points of the ordinary equator. 

The second region can be represented as a composition of two open sheets. One 
edge of each sheet is adjacent to the magnetic equator and the straight line passing 
through the center of the neutron star orthogonal to the rotational and magnetic axes. 
The other edge of the sheet extends to infinity, such that, for an arbitrary point lying on 
the force-free surface, 6" — )■ 7r/2 as r — t- oo. At large r, the sheet differs little from the 
6" = 71 /2 plane. However, the difference becomes appreciable at relatively small distances 
from the neutron star, and a cupola-like protrusion ending at the magnetic equator forms 
above the 6" = 7t/2 plane. The two sheets are joined smoothly at the points of the line 
9 = 6' = 6" = 71/2, forming a single sheet. The points of this line cannot be described 
using because this equation is degenerate at these points, in connection with our 
choice of coordinates on the surface of the angular coordinates {9,ip). We can directly 
verify that this line does indeed belong to the force-free surface using (jlj). 

As we can see from the figure, the dimensions of the closed arched parts of the force- 
free surface grow with 9m- When this angle reaches 7r/2, there is a break in the edges 
of the arches adjacent to the equator. These edges smoothly join with the edges of the 
open sheets that extend to infinity, forming an axially symmetrical, double-cupola figure 
specified by the equation = + 4cos^6''). It is easy to see that the electric field 
is also orthogonal to the magnetic field when 9 = 7i/2. Thus, the force-free surface of 
an orthogonal rotator is a combination of the figure described above and the equatorial 
plane. 

If the angle 9m is increased further, the part of the force-free surface that was arch- 
like for 9m < tt/2 is transformed into an open sheet when 9m > 7r/2, as a result of the 
edge breaking away from the equator and its transition to infinity. The open sheets, on 
the contrary, are transformed into arch-like parts of the surface, due to the attachment of 
their open edges to the equator. Note that we are speaking here of the breaking away of 
edges with some tentativeness, since it is more correct to speak of a reconnection of parts 
of the force-free surface as the angle 9m passes through n/2. Indeed, if 9m becomes equal 
to 7r/2, the closed arch-like parts become non-smooth while remaining continuous, and 
each of the arches is represented as a joining of the bulging part and a fiat part belonging 
to the equatorial plane 9 = 7i/2 and lying inside the double-cupola structure noted above. 
This is also true of the open sheets, with the exception that their fiat parts lie outside this 
structure. In the transition through 9m = 7r/2, there is a sort of reclosing of the parts of 
the force-free surface: the fiat part of the open sheet becomes reconnected to the bulging 
part of the arch, and the fiat part of the arch with the bulging part of the open sheet. 
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Taking this into account, we further consider everywhere only angles 6m lying be- 
tween and 7r/2. This does not limit the generality of our discussion, since a rotator 
with the angle x > 7r/2 between the vectors m and f2 is equivalent to a rotator with 
^m = 7r — X<7r/2, with Q replaced by — Q in all formulas. 



3. MOTION OF CHARGED PARTICLES 

We will investigate the motions of particles in a vacuum magnetosphere using the classical 
Dirac-Lorentz equation 



3c3 



X + —X X Xk 



c 



+ F\ (6) 



where = (ct, r-^)-^ is a contravariant four-vector containing the time t and coordinates 
r = {x,y,z)'^ of the particle in the laboratory frame (T denotes transposition), mg the 
particle's mass, e the particle's charge, and c the speed of light. A dot above a four- vector 
denotes differentiation with respect to the proper time r of the particle, i.e., the time in 
a frame comoving with the particle. The differentials dt of time in the laboratory frame 
and dr of the proper time of the particle are related as dr — dt/^, where 7 is the Lorentz 
factor of the particle. The four-force acting on the particle is given by 

c 

where Fik = dAk/dx"^ — dAi/dx^ is the electromagnetic- field tensor and the definition of 
the four-potential = {A'^ , A^ , A'^ , A^)'^ = (0, A"^)-^ contains the standard scalar and 
vector A electromagnetic potentials, with the electric and magnetic fields given by the 
formulas 

1 dA 

E = -V0 —, B = VxA. 

c ot 

The indices i and k take on the values 0, 1, 2, 3, with repeating indices everywhere denoting 
summation. The transition from contravariant to covariant components and vice versa is 
carried out using the metric tensor gik — g'^^ — diag(l, —1, —1, —1). 

We must obtain the equation of motion of the particle in the laboratory frame. This 
is straightforward using the following expressions for the derivatives of the four- vector 
with respect to the proper time of the particle: 

■ i dx' 

d'y dx"^ 9 d'^x^ ^rr\ 

^ (d2_Ydx^ od^dx^ 3^2^^ , .d^ 
^\dtJ dt dt'' dt dt dt'' dt^' 

Let us make a transition to a set of dimensionless variables. We will measure the 
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strengths of the electric and magnetic fields in units of the so-called critical field Bcr = 
mlc^/eh ~ 4.4 x 10^'^ G, the particle velocity in units of the speed of light c, the particle 
charge in units of the positron charge e, the particle mass in units of the electron mass mg, 
the particle energy in units of the electron rest energy mgC^, all distances in units of the 
Compton wavelength of the electron X = h/ rrieC ~ 3.86 x 10^^^ cm, and all times in units of 
X/c. Note that, in these units, 1000 km ^ 2.6 x 10^*^ and 1 s 7.8 x 10^°. The convenience 
of using these units is due to two factors: the presence of strong electromagnetic fields in 
the magnetosphere, which makes it reasonable to measure the field strengths in terms of 
the critical field, and the ultra-relativistic motion of the particles, which makes it natural 
to characterize the particle energy in terms of its Lorentz factor. 

Using ([7]) after introduction of the above dimensionless variables and separation of 
the scalar and vector components of the four- vector x*, the Dirac-Lorentz equation (jSj) 
reduces to the system of equations 



d'j 
Itt 


2 

= 3«7 


df^ ' \dtj _ 


± V ■ E, 


dv 
^~dt 


2 

= 3«7 


^ dj dv d^w 
dt dt dt"^ 


± 


E - v(v • E) + V X B 



where v is the particle velocity and a = /he ~ 1/137 — is the fine-structure constant. 
When considering the motion of positrons, we should take a "+" sign in the system ([8]) and 
0, and should take a "— " for the motion of electrons. For convenience, we will consider 
the "+" case, and comment when necessary on the changes that result if electrons rather 
than positrons are considered. 

Equation ^ represents the conservation of energy, and the equation of motion 
of the particle. This is easy to see if we formally specify the fine-structure constant to be 
a = in dH]) and ([9]). These equations then take on their standard form, corresponding to 
neglecting the inverse infiuence of the field of the moving charged particle on the particle 
itself. The first of these equations indicates that the energy acquired by the particle per 
unit time is equal to the work done by the electric field on the particle. The second of 
these equations, ([9]), simply refiects the fact that the change in the relativistic momentum 
of the particle, p = 7V, per unit time is due to the action of the total Lorentz force, 
E + V X B. 

One might ask why we are using the classical Dirac-Lorentz equation to describe 
particles in the electromagnetic field of the vacuum magnetosphere of a neutron star, and 
not the usual equations of motion. The reason is that the electric field is so strong that 
electrons and positrons are virtually instantaneously accelerated to relativistic energies, 
and there is intense emission of so-called curvature radiation when these particles move 
along the curved magnetic-field lines, due to their very high Lorentz factors, leading to 
energy losses by the particles. This radiative friction force is taken into account by the 
Dirac-Lorentz equation. The use of the usual equations of motion would lead to ap- 
preciable over-estimation of the particle energies, and therefore also of the energies of 
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the photons emitted by the particles. In studies of amphfication processes in the elec- 
tron-positron plasma in the vacuum magnetosphere, this would lead to results that bear 
no relation to reality. The presence of radiative friction must be taken into account when 
numerically computing the trajectories of charged particles in the inner magnetosphere 
of a neutron star It is interesting that, in vacuum fields, this is also important at 



distances of the order of the light cylinder [27|, |28|, when considering the behavior of 
charged particles in a wave field and a constant electric field 29 1. 

Let us consider this question in more detail. A charged particle created in the 
magnetosphere, be it an electron or positron, will experience an electrical force and be 
accelerated. Let us estimate the time required for the particle's motion to become rela- 
tivistic. For this, it is sufficient to use the equation of motion of the particle projected 
onto the direction of the magnetic field, dp\\/dt = E^\. Since the electric field does not 
change significantly over the time the particle is accelerated, we can immediately obtain 
the characteristic time Trei for the particle to make a transition to the relativistic regime, 
when the longitudinal momentum becomes close to py ~ 1: 

Trel ^ (10) 

We see from the general form of ([T]) and (j2]) how the electric field E is related to the 
magnetic field B: 

where the radius of the light cylinder in ordinary, dimensional units is Rl = c/Q, and 
takes on values for typical pulsar periods of P ~ 0.1—1 s of Rl ~ 10^—10^ km (we have 
used the fact that the wave number is A; = 1/Rl). We can find the electric field at the 
stellar surface E^urf by setting r = R: 

R 

Esurf — PsMr/' (H) 

For a typical surface magnetic field Bg^rf ~ 0.01—0.1 and a ratio of the radii of 
the neutron star and the light cylinder R/Rl ~ 10~^— 10~'^, the surface electric field is 
Esurf ~ 10~^— 10~^. Accordingly, the time for the transition to the relativistic regime 
is Tj-e; ~ 10^—10^. This means that the particle will reach near-light speeds after a 
time of the order of 10~^^ — 10~^^ s, after which its motion can be taken to be ultra- 
relativistic. Obviously, the particle will traverse a distance of no more than 10^—10^ 
Compton wavelengths during this time. We have obtained an upper limit — the particle 
velocity has not reached the speed of light c in the initial stage of the acceleration in the 
time interval considered. This justifies our assumption that the electric field does not 
change significantly during the particle's acceleration time, since the acceleration time is 
small compared to the period of rotation of the neutron star, and the distance over which 
the acceleration occurs is small compared to the distance over which the electric field 
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changes appreciably; in our case, this is a distance of the order of the radius of the star, 
~ 10 km. 

Further, the particle will continue to be accelerated in an ultra-relativistic regime. 
As is well known, in general, the particle will move along the curved magnetic-field lines. 
Acquiring more and more energy, the particle begins to emit curvature radiation, whose 
characteristic energy Scurv and total intensity Wcurv are equal to 



37^ U/ -I l' 

o ; curv o 

2 p 3 p 



W,urv = Tja^i' (12) 



where p is the radius of curvature of the particle trajectory, the energy Scurv is measured 
in units of mgC^ (like the particle energy), and the intensity Wcurv is measured in units 
of rrieC^/X. Since the intensity of curvature radiation grows with the Lorentz factor 7, 
the particle will eventually not undergo further acceleration after it has reached some 
maximum Lorentz factor 70, since all the energy it acquires from the electric field will be 
lost to curvature radiation. To determine 70, we must find a stationary solution to OH]), 
having substituted in this equation d'j/dt = d?'-)/dt^ = 0: 

%j(5)^ = v.E. (13) 



3 "\dt/ 

Equation f|T3|) shows that, in the stationary state, all the work done by the electric field on 
the particle is completely transformed into the energy of curvature radiation (fT2l) . In the 
case considered, the particle's velocity vector has unit length, f = 1, so that dv/dt = n/p, 
where n is the principle normal vector to the particle trajectory and p is the radius of 
curvature of the trajectory. The maximum Lorentz factor of the particle takes the form 

/ X 1/4 

7o=(^i?||P^) . (14) 



For a characteristic longitudinal electric field E\\ ~ 10~^ and radius of curvature of the 
particle trajectory p R c:^ 2.6 x 10^^, the maximum Lorentz factor of the particle is of 
the order of 70 ~ 6 x 10^. Here, we have taken as a characteristic radius of curvature of 
the particle trajectory the radius R of the neutron star, having in mind that the particle 
moves along trajectories close to the magnetic-field lines, and the radius of curvature of 
these lines is of the order of R near the magnetic equator. 

Let us estimate the time for a particle to acquire its maximum Lorentz factor 70 — the 
time for the full acceleration of the particle and its transition to a quasi- stationary motion 
regime determined by the equilibrium condition (fT3|) . We use here the conservation of 
energy of the particle in its simplest form, d'j/dt = E\\, without including the radiative- 
friction force. The desired time is then 
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We can see from ffTOjl and ffTSjl that Tgt = To^rez; i-e., the time for the acquisition of the 
maximum Lorentz factor is a factor of 70 greater than the time for the particle to achieve 
near-hght speed. Thus, the particle acquires its maximum energy over a time of the order 
of Tst ~ 10^2 (in dimensional units, 10~^ s), after which the work by the electric field 
on the particle per unit time is equal to the total intensity of the curvature radiation. 
During this time, the particle traverses a distance of the order of 10^^ (in dimensional 
units, several tens of centimeters), appreciably less than R; therefore, the assumption 
that the electric field does not change significantly over the total acceleration time Tst is 
valid. This means that we can take the acceleration of the electrons and positrons to be 
virtually instantaneous, and to occur at the point where the electron-positron pair was 
created. 

However, both the electric field strength and the radius of curvature of the trajectory 
change during the motion of the particle, so that the Lorentz factor 70 varies with time. 
This means that the particle will adjust its motion, acquiring energy due to work by the 
electric field if 70 grows along its trajectory or losing energy to curvature radiation if 70 
decreases along its trajectory. It is now important to determine the rate at which this 
adjustment of the particle's energy occurs. If this rate is appreciably higher than the rate 
of variation of 70 along the trajectory, we can take the Lorentz factor to be determined 
by the coordinates above the point where the particle is located; at a fixed moment in 
time, the value of 70 depends only on the coordinates of the point considered, not on the 
velocity of the particle. 

Let us find the variation of the Lorentz factor 7 of the particle as it approaches 
the steady-state value 79. Here, we use ([H]), having represented the Lorentz factor of 
the particle as the sum of the steady-state value 70 and some deviations from this value 
57. The deviations ^7 can be taken to be small, ^7 ^ 70. This is true because we are 
considering how the Lorentz factor of the particle changes during its motion as a result 
of the smooth variation of 70 due to the variation of the longitudinal electric field E\\ and 
the radius of curvature of the trajectory p, with the particle initially being in a stationary 
state and possessing an energy 79. Using the smallness of the deviations ^7 to linearize 
([8]), we obtain 

The characteristic equation — 3(j/2«7o — ^ll/ = for (fT6|l has the solution 

3 87!^ 

^1 = 7. , ^2 = --a^- 17 

2q;7o 3 

When computing these equations, we used the condition |cr2/cri| ^ 1. It follows from 
()14l) and f[T7|) that the condition that the ratio of the numbers themselves be small is 
equivalent to the condition 

\aE\\ < L (18) 
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The condition f|T8l) is always satisfied, and the indicated quantity of the order of 10~^ 
for the characteristic longitudinal electric field, E\\ ~ 10~^. Note that the inequality f[T5|) 
does not represent a constraint, and is satisfied even for the critical electric field E\\ = 1 
(it is well known that the electric-field strength cannot exceed this value due to the direct 
creation of electron-positron pairs from the vacuum in the presence of such strong fields) . 

The general solution of (I16p takes the form 

5^ = Cie'i* + 026"^'. (19) 

As we can see from f[T7|) . the eigenvalue cxi is positive. Formally, this would mean that a 
particle given some additional energy above 70 would begin to accelerate further, acquiring 
more and more energy. Clearly, this cannot happen physically. This is the so-called 
paradox of self-acceleration of the particle, well known in the theory of the Dirac-Lorentz 
equation [30j. It arises because the Lorentz-Dirac equation ([6]) contains a third derivative 
of the four-vector with respect to the proper time of the particle r. This means that 
it is not sufficient to know the initial coordinates and velocity of the particle in order to 
describe its motion, so that boundary conditions must be specified for this equation. These 
boundary conditions must be chosen to eliminate the self-accelerating solution. Apart 
from specifying the initial coordinates and velocity, we must require that the acceleration 
of the particle becomes zero after all external forces cease to act on it. In our case, 
this corresponds to having the acceleration of the particle vanish as t — ?■ 00. For the 
general expression ( 1T9|) . this condition will be satisfied if and only if Ci = 0, so that the 
self-accelerating solution disappears. 

Taking this into account, we can immediately conclude that, when the Lorentz factor 
of the particle deviates by an amount (57i from the stationary value 70, the Lorentz factor 
will approach 70 according to the exponential law 

57 = 57. e-*/"« 



with the decay constant 

^0 = f 4- (20) 

8a 7o 

For characteristic radii of curvature of the trajectory p ~ 10^'' and Lorentz factors 
7o ~ 10^— 10*^, the decay time is tq ~ 10^^ — 10^^ (in dimensional units, 10^^'^— 10^^ s). 
Over a time Tq, the particle travels a distance /q = '^"0 of the order of 10^^—10^^ Compton 
wavelengths (i.e., from centimeters to several tens of meters). Now recall that the charac- 
teristic distances over which the electric field and radius of curvature of the magnetic-field 
lines — and therefore also the stationary Lorentz factor 70 — vary are of the order of R. Be- 
cause /o ^ particles have time to adjust to variations in 70 during their motion in the 
magnetosphere. Thus, we can take the Lorentz factor of a particle to be determined fully 
by its coordinates. 
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This assertion requires some refinement. We have aheady discussed the fact that 
there exists a force-free surface (jS]) in the magnetosphere. It is clear that charged particles 
will move toward this force-free surface. However, the longitudinal electric field will 
gradually decrease as a particle approaches this surface. Simultaneously, the time for 
the adjustment of tq will increase. Eventually, the longitudinal electric field becomes 
sufficiently weak that the reaction time of the particle is too great to satisfy the condition 
of quasi- stationary motion; i.e., to support a balance between the work done by the 
electric field and the intensity of the curvature radiation. This is associated with the fact 
that the particle Lorentz factor decreases as the field weakens, together with the energy 
and intensity of curvature photons. As a consequence, the characteristic time for the 
particle energy losses near the force-free surface begins to exceed the characteristic time 
for variation of the electric field. We will consider this question in detail in our analysis 
of the capture of particles by the force-free surface. 

Thus far, we have discussed the energetics of particles using only (|8]). Let us now 
consider in more detail the second equation (|9]), and investigate how this equation affects 
a particle's trajectory. We will estimate the magnitudes of terms arising due to our al- 
lowance for the self-interaction of the charged particles, assuming that a time Tst ( |T5|) has 
passed after their creation, and that the particles have by this time already fully accel- 
erated and made a transition into a quasi- stationary regime determined by the condition 
(fT3!) . Since the particle energy is determined by the Lorentz factor 70 (IT^ . which varies 
over distances of the order of the radius of the neutron star R, d'-f/dt ~ 7o/-R- The accel- 
eration of the particle is of the order of \d\/dt\ ~ 1/-R, since the radius of curvature of the 
particle's trajectory is close to the radius of curvature of the magnetic-field lines, which, 
in turn, is of the order of R, if we consider a region near the magnetic equator of the star. 
According to the same reasoning, the second time derivative of the particle velocity will 
be Id'^v/dt'^l ~ 1/R^. Thus, the first term in ([H]) containing the square brackets is of the 
order of a'-^^/R^. This is small compared to the second term, since 

«^ « (21) 

as we can see by taking typical values of the Lorentz factor 70 ~ 10^, the stellar radius 
R ~ 10^'', and the electric field E ~ 10~^. However, the fulfilment of the condition (12T|) 
is not a sufficient basis to neglect terms allowing for the self-interaction. The term on the 
left-hand side of OH]) is to order of magnitude 70/-R ~ 10~^, and is also appreciably smaller 
than E. To justify neglecting terms allowing for the infiuence of the fields created by the 
charged particles on their own motion, these terms must be small compared to 70/-R: 

We see that this condition will always be satisfied, so that the motion of the particles in 
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the vacuum magnetosphere of the neutron star will be described by the equation 

7^ = E- v(v ■ E) + V X B. (22) 

This equation has been written for positrons; the corresponding equation for electrons 
will differ from (122|) only in the presence of a minus sign before the entire right-hand 
side. Thus, Eq. f l22|) describes the proper motion of the particle, while its energetics are 
determined by ([H]), which reduces to f[T^ in the quasi- stationary case. 

Multiplying ( 122|) vectorially by B, we obtain an iterative formula for — the velocity 
component orthogonal to the magnetic field: 



1 



E- v(v-E) -7— 

^ ^ ' dt 



X B. (23) 



Generally speaking, charged particles in the strong magnetic field of a neutron star 
will have virtually no momentum orthogonal to this field: due to synchrotron cooling, 
the particles immediately go to the zeroth Landau level. This means that, in a first 
approximation, we can assume that the particle velocity is directed along the magnetic 
field, V = b, where b = B/i? is a unit tangent vector to the local magnetic-field line. 
In a slightly more refined approximation, we must take into account the presence of a 
drift-velocity component orthogonal to the magnetic field. We will search for solutions in 
the form v = b + v_|_. Note that the quantity v_|_ is first order in R/Rl. After substituting 
the relation v = b into (!23l) . the total velocity of the particle v can be written 

V = b + Ve + Vc, 

where Vg is the electric drift and Vc the centrifugal drift velocity. These velocities are 
defined by the formulas 

ExB 

= (24) 

V, = 1 b X ^. (25) 
B dt ^ ^ 

It is easy to write the total time derivative of the tangent vector to the local magnetic-field 
line as 

dh dh , 

— = — + V - V b. 

dt dt ^ ' 

The first term on the right-hand side is to order of magnitude \dh/dt\ ~ f2, or 
equivalently ~ 1/P, and, for typical pulsar periods P ~ 0.1—1 s takes values 10~^^ — 10^^°. 
Due to the condition v = b, the second term is (b ■ V)b = n/p, where n is the principle 
normal vector to the magnetic-field line and p the radius of curvature of the field line. For 
typical radii of curvature p ~ 10^^, the first term can be neglected compared to the second, 
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so that the velocity f l25|) mainly describes the centrifugal drift, which can be written 



V, 



c 




(26) 



where e = b x n is the binormal vector. Note the difference in the directions of the 
centrifugal drift velocities for the electrons and positrons: formula f l26p is written for 
positrons, and for electrons, we must add a negative sign before the right-hand side. The 
electric drift velocity Ve for both electrons and positrons is given by ( 12^ . and the direction 
of this velocity does not depend on the sign of the charged particle. 

Let us estimate the drift velocities Vg and Vc to order of magnitude. The electric 
drift velocity is equal to the ratio of the electric and magnetic fields, Ve ^ E/B, and is of 
order R/Rl ~ 10"^ [see ([II])]. For 7 ~ 10^ p ~ 10^^ and B ~ 0.01-0.1, the centrifugal 
drift velocity is Vc ~ 10^^—10"'^. Thus, the centrifugal drift velocity is the next order 
of smallness compared to the electric drift velocity, and we will accordingly neglect the 
centrifugal drift. We will assume that 



with accuracy to within terms o{kr), i.e., to within quantities of quadratic and higher 
order in kr ~ R/R^. Note that, for fl27|) . the equality f = 1 is satisfied in this case with 
accuracy to first order in R/R^, and a discrepancy appears only in second order. If it 
were necessary for us to take into account quantities of the order of ~ 10^*^, we would 
have to calculate the following terms in the expansion of order (R/Rl)^ by iterating f l23|l . 
This would lead to corrections to the longitudinal velocity of the particle, and it would 
no longer be possible to assume vy ^ b. The satisfaction of ( 127|) and the smallness of 
the drift for the orthogonal motion compared to the longitudinal motion enabled us to 
assume that the particles essentially move along the magnetic-field lines. 

The electromagnetic field described by (II]) and ([2]) is periodic in time. This leads 
to a certain inconvenience when considering the motion of a particle in the laboratory 
frame. If the particle travels some distance dl between points 1 and 2 over a time dt, 
the variation of the electromagnetic field at the final point 2 is determined not only by 
the variation in the coordinates, but also by the fact that the field itself at the point 
2 varies over the time dt. Let {t, r, 6, (p) be the spherical coordinates in the laboratory 
frame. We transform to the new coordinates (t', r', -(9', (p') using the relations t' = t, r' = r, 
i)' = 6, if' = p — Qt. This leads to the transformation of the partial derivatives d/dt = 



d/dt' -nd/dip', d/dr = d/dr', d/dO = d/d^', d/dip = d/dp'. We will call {f , r' , ^' , <p') 



the coordinates in the rotating frame. In these coordinates, the electromagnetic field 
depends only on (r', i)', p'), not on t'. Further, we will omit the primes for variables when 
it is clear that we are considering quantities in the rotating frame. The transformation of 
the velocity and acceleration in the transition to the rotating system have the standard 



V = b + Ve 



(27) 
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form 



V = v' + Vtr, 




(28) 



■r, 



where v' is the relative vefocity and vtr = ri x r the translational vefocity. Formulas (128|1 



Here, E and B are the electric and magnetic fields, which are determined by ([2]) and ([T]), 
but depend on the coordinates (r', i)', cp') after the described change of variables. 

When obtaining (129|) . we neglected the Coriolis, 2 x v', and translational, ft x vtr, 
accelerations compared to the relative acceleration dV/dt'. This is possible because the 
relative acceleration is primarily axipetal, and due to the motion of the particles along the 
curved magnetic-field lines. Therefore, \d\'/dt'\ ~ l/i?, which is to order of magnitude 
10^^^. To order of magnitude, the Coriolis acceleration is |2r2 x v'| ~ which is 

10^^^. The translational acceleration is still smaller, because the translational speed Vtr 
is clearly less than t;' ~ 1 in the considered regions in the magnetosphere, r < lOR Rl, 
where the efficient single-photon creation of pairs is possible. Moreover, only when r ^ R^ 
are ([1]) and ([2]) for the electromagnetic field valid. In general, in this case, the ratio of the 
translational to the Coriolis acceleration is equal to the ratio of the Coriolis to the relative 
acceleration, and is R/Rl ~ 10~^. Thus, the neglected term in the left-hand side of ( 12^ . 
7f2 X (2v' -|- Vtr), is also small compared to the electric-field strength, E ~ 10~^; for 
7 ~ 10^, it is to order of magnitude "j/Rl ~ 10~^^. In the right-hand side of (l29i) are left 
only terms whose magnitudes are comparable to E, while all substantially smaller terms 
are neglected. The term v' x B is retained, because, although f' ~ 1, the velocity vector 
itself v' is nearly parallel to the magnetic field, and its orthogonal component v^, which 
is the only one to contribute to the vector product, is to order of magnitude Ve ~ E/B, 
so that v' X B is also of the order of E. The Lorentz factor 7 in fl2^ is measured in the 
laboratory frame. 

Further, we will investigate the motion of the particles in the rotating frame. We are 
interested in whether there exist regions in the magnetosphere where the accumulation of 
primary plasma is possible. We will find equilibrium positions — points where a charged 
particle can remain for an indefinitely long time. The coordinates of the equilibrium 
positions are determined by the conditions that the velocity and acceleration be equal to 
zero in the rotating frame: v' = dV/dt' = 0. Substituting this condition into (12^ leads 
to the equation 



were we have introduced the effective electric field 'E^^^ = E + v^^ x B. The components of 



are valid when Cl Q^. This is equivalent to the condition that P ^ 1, which is always 
satisfied, since, to order of magnitude, P ~ 10"^^. Substituting fl28l) into fl22l) yields for 
the equation of motion of the particle in the rotating frame 



7^ = E + X B - v'(v' • E) + v' X B. 



(29) 



E'^f = 



(30) 
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this field, which can straightforwardly be found from ([T]), ([2]) together with the expression 
for \tri have the following form in spherical coordinates 



Et^f = -kr 



— COS OBr + (1 T ) sin OBg 

j - sin OB^ + cos 6Bq 
Ef = kr{l--^YoseB^, 



(31) 



where 5^, Bq and B^p are determined by ([T]). Substituting the condition (130|) into ([3T 
yields the set of equilibrium points 



vr 

r = R, ^ ~ 2 (equator), 
vr 

r = R, 9' = — (magnetic equator). 



(32) 



ri 3 + cos Or, 



1 — cos 9 m 



(open sheets) 



2 " 7 ' V 2 



(33) 



3 - cos ^„ 



i?2 1 + cos^^ 



(folds). 



(34) 



We can see that all the equilibrium points are located on the force-free surface ([S]) and 
can be divided into two groups: non-isolated and isolated. The nonisolated equilibrium 
points are all on the equator and magnetic equator [see f l32|) ]. Two isolated equilibrium 
points are added to these regions: two whose coordinates are given by (ES]), located on 
open sheets of the force-free surface, and two others given by on the folding parts of 
the surface adjacent to its boundary with the equator and magnetic equator. All four of 
these points lie in the plane passing through the rotational axis and magnetic axis. Only 
the points (!33|) lie in the cross section of the bisectrix angle between fi and m with the 
force-free surface, while the points f lM|) lie on the cross section of the line orthogonal to 
this bisectrix lying in the plane of the vectors O and m with the same force-free surface. 
The radial coordinates r_ of the points lying on cupolas are always smaller than the radial 
coordinates r+ of the points lying on open sheets of the force-free surface. The coordinates 
r+ and r_ are equal, r+ = r_ = -\/3_R, only if 9ra = 7r/2. If 9m — ^ 0, then r_ ^ R and 
r+ — oo; the equality r+/R = 2\/2/ 9m is asymptotically obeyed for small values of 9m- 

This raises the question of the stability of the equilibrium positions we have found. 
Moreover, the character of the trajectories of the particle motion is not clear, both near 
the force-free surface in general and near the equilibrium points in particular. We will 
provide an answer to this question after a detailed study of the capture of particles by 
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the force-free surface. 



4. OSCILLATIONS OF CHARGED PARTICLES NEAR THE 

FORCE-FREE SURFACE 

Let us consider the motion of charged particles near the force-free surface, applying some 
qualitative reasoning. We can see immediately that the motion will be oscillatory, because 
the longitudinal electric field changes its sign in the transition through the force-free 
surface. Here, we will assume that the signs of the charged particles are such that the 
electric force on them is directed toward the force-free surface. Further, some regular 
motion along the force-free surface will be superposed on the oscillatory motion of a 
charged particle. If we imagine carrying out an average over the rapid oscillations, thereby 
distinguishing the leading center, or the local equilibrium position about which the particle 
oscillates, this center can only be located on the force-free surface. Otherwise, on average, 
a non-zero electric field would act on the particle, causing it to return to the surface. 

We will study the particle motion quantitatively in the rotating frame using fl2^ . We 
choose some point fq on the force-free surface and expand the electric field E, magnetic 
field B, and effective electric field E^'^^ introduced above about this point as follows: 

E = Eo + (x' ■ V)Eo, 

B = Bo + (x' ■ V)Bo, (35) 
E'" = E'J^ + (x' ■ V)E^^^, 

where Eg, Bq, and Eq^'^ are the field strengths at the point tq, and x' = r — fq is the 
distance from the point r for which we are interested in the field strengths to the point 
fq. The expansion ( 135|) is valid because we are considering the motion of a particle in the 
immediate vicinity of tq, at distances appreciably smaller than the characteristic distances 
for variations of the fields, so that x' <^ R. Before searching for the full solution of f l29|l . 
we will first find some partial solution, describing the motion of a particle with constant 
velocity Vq = const. This solution will then satisfy the system of equations 

dt' ~ °' (36) 
= E^^^ + (x'o • V)E^^^ + v; X Bo + v; X (x[, ■ V)Bo. 

When obtaining the second equation, we used the fact that dw^/dt' = 0. We also required 
that Vq <^ 1, and neglected the term Vq(vq-E), which is small compared to E, and therefore 
also compared to E'^-^'^ . It follows from the first equation of fl36|) that Xq = v'^t (at the 
initial time, the particle is located at the point ro on the force- free surface). Substituting 
Xq into the second equation yields 

E^^^ + v[, X Bo + t 
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v)e; 







+ Vq X 



V)Bc 



For the motion of the particle to occur with the same constant velocity at the next moment 
in time, two conditions must be satisfied: 



E^^^ + X Bo = 0, 
[v[,-V)E^^^ + v[,x K-V)Bo = 0. 



(37) 



The first equation of 0371) unambiguously determines the velocity component orthog- 
onal to the magnetic field 

^ (38) 



S2 



with the longitudinal component v',, remaining arbitrary: 



(39) 



The form of fl38l) for the velocity coincides with that of expression flMl) for the electric 
drift velocity in the laboratory frame, with E replaced by E^-'^'^. It is easy to verify 
using 0381) that the drift velocity in the rotating frame is indeed given by the sum of 
the electric drift velocity in the laboratory frame and the "translational" drift, equal to 
b(vjr -b) — = — _L, which is simply the negative of the component of the translational 
velocity orthogonal to the magnetic field. This is true because we are considering the 
motion in the rotating frame; the presence of the translational velocity Vf^ = x r provides 
an additional contribution to the drift, so that, in the laboratory frame, Vg = + vtr±- 
We introduce the velocity components v^: 



1 



■'±r 



fir sin 9 brb^- I 1 



fir 



cos 9 br 



1 ) sin 9 ba 



-fir 



sin 6* ^-6^ + 6gj + cos 9 brbg 



where br, bg, and b^ are components of the unit vector b. 

The longitudinal component tiy of the velocity Vq is unambiguously determined by 
the second equation of fl371) . Taking the scalar product of this equation and the vector 
Bo yields 

" (Vn 



Using the fact that Eq- 



Bo 

eff 



,o-V)E^^-^ + (BoX v[,)-(v[,-V)Bo = 0. 
-Vq X Bo, we immediately obtain 



'0 



VfEo ■ Br 



0. 



(40) 



Here, we have also used the fact that Eq-^-^ 



Bo = Eq • Bq. However, since the gradient 
V(Eo ■ Bo) is directed normal to the force-free surface, E ■ B = 0, the velocity Vq lies in 
the tangent plane passing through the point tq of the force-free surface. 
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Thus, the partial solution describing the motion of a particle with constant velocity 
exists and is given by f l5^ . The longitudinal velocity component is obtained from the 
condition (Hn|l . and has the form 

^, ^ ■ V(Eo ■ Bo) 



b ■ V(Eo ■ Bo) 



We can see that the velocity Vq lies in the tangent plane, so that the particle cannot 
leave the force-free surface, and its trajectory lies entirely on this surface. The velocity 
will vary from point to point, and is determined by the same expression f l39p . but with 
the fields at the point where the particle is located at the given time used in place of Eq-'^''^ 
and Bo. 

However, the solution found does not exhaust all classes of motion of the particle. 
This is clear, if for no other reason than because a solution with velocity Vq describes 
adiabatic motion of the particle along the force-free surface, with no oscillations of the 
particle whatsoever. We will now find the full solution of (12^ near the force-free surface. 
We seek a solution in the form 

x' = x[, + x;, v' = v[, + v;, (41) 

where Xq and Vq are the adiabatic solution found above. No special restrictions are applied 
to the quantities x'^ and v'^^ a priori, apart from the requirement that x'^^ be small compared 
to the characteristic scale for variation of the electromagnetic field, R. In particular, the 
velocity v'^ can be close to the speed of light. Substituting ( ]¥T]) into (129]) and using (15^ 
and f l35|) leads to the equation 



7^ = {< ■ V)E^^^ + v[, X (x; ■ V)Bo + v; X B - ■ E), (42) 

where E and B are determined by (!35|) and the coordinate value x' by fj4T]) . When 
deriving fH2|) . we used the condition v'q <^ 1. The term v'^(v']^ ■ E) is kept because, 
generally speaking, the velocity v'^^ can be of the order of unity, so that this term can be 
comparable to E. 

Let us first consider the non-relativistic case {v[ 1). It immediately follows from 
(USD that 

^ = (x; ■ V)E^^^ + V, X (x; • V)Bo. (43) 

Here, we have used the fact that v'^^ ^ v[h, i.e., the velocity is directed essentially along the 
magnetic field, so that the term x B vanishes. It is obvious that the linear differential 
equation ( 143|) describes oscillations. The left-hand side contains the second time derivative 
of the vector x'^^, and the form of the right-hand side is linear in the components of the 
vector x'^, and can be represented as the product of a 3 x 3 matrix and the vector x^ 
However, we can convince ourselves of the oscillatory character of the solution without 
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searching for the eigenvalues of this matrix and corresponding eigenvectors. Let us find 
the frequency for the non-relativistic oscillations. We take the scalar product of fH51) 
with the vector Bq. Carrying out manipulations on the right-hand side of the resulting 
equation analogous to those described above for the derivation of fHOl) straightforwardly 
yields for the square of the frequency of the non-relativistic oscillations 



= _^L^(^^, (44) 
Bo 

In this b and v'^^ ^ v[h, and the coordinate x[ satisfies the oscillation 

equation (Px'i/dt^ = —u'^x[. When deriving OS]), we also assumed implicitly that the 
instantaneous direction of the vector b at the point where the particle is located coincides 
with the direction of the vector bo at the point fq about which we expanded the fields, 
i.e., bo ~ b. 

Let us find criteria for applicability of the non-relativistic approximation. As we can 
see from fH4|) . the characteristic oscillation frequency is a; ~ VOB, or equivalently. 



Here, we have used the estimate (fTTj) for the electric field at the surface of the neutron 
star. For characteristic light-cylinder radii Rl ~ 10^° (for pulsars with P ~ 1 s) and 
magnetic fields B ~ 0.01—0.1 the non-relativistic oscillation frequency is of order u ~ 
10^^^— 10^^*^ (which corresponds to u = ujl-n ~ 1 — 10 GHz in dimensional units). In the 
non-relativistic case, the v\ <^ 1, and the oscillating particle travels a distance that is 
clearly less than l/u, during one period. The maximum amplitude of the non-relativistic 
oscillations is ^ 

Inro — • (^^) 
OJ 

The oscillatory motion will be non-relativistic if x'l <^ Inro- For the frequencies u found 
above, we have Inro ~ 10^°— 10^^ (0.1—1 cm in dimensional units). This means that the 
equality bo = b used when deriving fl44p is precise by virtue of the extreme smallness of 
Inro compared to R. We can see that oscillations of a charged particle about the force- 
free surface with an amplitude exceeding 1 cm are clearly relativistic. Thus, our primary 
interest is in considering the relativistic case of oscillations, which is realized in nature. 

To obtain an equation describing the particle oscillations without placing any special 
constraints on the velocity, we must use fH2l) . setting in this equation || B: 

7^ = (x; ■ V)E^^^ + v[, X (x; ■ V)Bo - v;(v; • E). (4?) 
Taking the scalar product of fjTTj) with Bo yields 

7^ = -^M-(.0^i?||, (48) 
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where u is the frequency of the non-relativistic oscillations given by fH4|) . We can see 
immediately that the term {v[YE\\ becomes negligible and 7^1 when v'^ ^ 1, so that 
(HSj) is transformed into the usual equation for non-relativistic oscillations with frequency 
u. Now let the velocity v'^ not be small. We must then use the condition x'^ <^ R and 
write the longitudinal electric field in the form E\\ = —ij^x\. After using the relation 



7 1/a/1 — (fi)^, we quickly obtain the equation for relativistic oscillations of a charged 
particle about the force-free surface: 

We will first make a number of qualitative comments concerning the ultra-relativistic 
regime for the particle's motion, when its velocity is close to the velocity of light, v'^ = 1, 
during a large fraction of the oscillation period. Let us find the regions where the particle's 
motion is non-relativistic. It is clear that these regions are located near the turning points 
of the motion, where the particle's velocity vanishes and its positional coordinate reaches 
the amplitude A. The width of the regions of non-relativistic motion near the turning 
points is then determined by the so-called acceleration length la, equal to the distance over 
which a particle that is initially at a turning point and has zero velocity acquires near-light 
speed. Another way of expressing this is to say that, after traveling a distance equal to the 
acceleration length the relativistic momentum of the particle p becomes comparable to 
unity. An upper limit for la is given by the amplitude Inro for non-relativistic oscillations 
of the particle [see Eq. ( H6|) ]. Indeed, if a particle that is initially at rest deviates by the 
distance Inro from the equilibrium position Xq and attains relativistic energies when it 
passes through the force-free surface, a particle located at a turning point will clearly be 
accelerated over a shorter time. This is true because the electric field at the turning point 
is non-zero, and exceeds the nearly-zero field in the region of non-relativistic oscillations. 
A lower limit for la is given by the characteristic time for the acceleration of particles to 
relativistic velocities in the magnetosphere r^ei (see ( !T0|) ). Thus, Trei ^ /a ^ Inro- It is 
easy to obtain a more accurate estimate of the acceleration length la near a turning point 
using the same reasoning as was used when deriving fllU]) . but setting E\\ ~ u'^A: 

A ■ 

The ultra-relativistic approximation is applicable if la ^ A. Note that the acceleration 
length la is smaller than Inro by a factor of A/ Inro, i-e., appreciably less than 1 cm. This 
obviates the need to consider the exact character of the motion of a charged particle near 
its turning points for oscillation amplitudes of the order of several Inro- In this case, the 
particle moves from one turning point to the other with a velocity v[ = 1, so that the 
oscillation period is T = 4A. The oscillations themselves are given approximately by 

Xi = — arcsm sm I — ) (50 j 
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and represent saw-like oscillations with amplitude A, with x'l = t when —A ^ t ^ A 
(further, we do not distinguish t and t', since t = t'). 

Nevertheless, ( H^ can be analyzed exactly without using the ultra-relativistic ap- 
proximation. This equation has the first integral 

c = . + ^. 

This first integral is simply equal to the maximum Lorentz factor 7max reached by the 
charged particle in its motion through the force-free surface; i.e., C = 7max when x'l = 0. 
The first integral C can be expressed in terms of the oscillation amplitude: 

C^l.^. (51) 

It can be shown that, after the following changes of variables. 



sm0o = sm 



(52) 



C-1 C-1 

Eq. (jl9j) can be written in the differential form 

'/2C 1 . 



dt = ± dlEU, k) - cos>o^(0, ft:)!- 

u sin 00 

Here, F{(f), k) and -E(0, k) are elliptical integrals of the first and second kind, respectively. 



defined as is done by Gradshteyn and Ryzhik 3l| : 

, . . da' 







1 — sin^ a' 

(53) 



E{^, k) = y a/i - ^2 sin^ a' da'. 



Introducing the function -R(0, k) = £'(0, k)— cos^0o-^(0) we obtain the relation between 
(j) and t 

t = ^^(M. (54) 

iij sin 00 

We can immediately exactly determine the oscillation frequency, noting that = cor- 
responds to the equilibrium position and = 7r/2 to a turning point. The distance from 
the equilibrium position x'^ = to the turning point x'^ = A is traversed by the particle 
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over one-quarter of its period, so that 



T = 4^m, (55) 

u sm (pQ 

where 'R{k) = E(k) — cos^0oK(k), and K(k) = F(7r/2, k) and E(ft) = E{tt/2, k) are total 
eUiptical integrals of the first and second kinds. 

Let us verify the asymptotic oscillation period fl55|) . If the oscillations are non- 
relativistic {A ^ Im-o), as we can see from (l5Ti) and (l52i) . C ~ 1, k ~ 0, sin0o — 1/^2, and 
K(0) = E(0) = tt/2, so that T = 27r/u; is the usual period for non-relativistic oscillations. 
If the oscillations are ultra-relativistic {A ^ Inro), then C ^ uj'^A'^/2, h ~ sin 00 — 1, 
and E(l) = 1, so that T = AA. We obtained this same result above based on qualitative 
reasoning. Note that, in spite of the logarithmic divergence of K(k) as k — )■ 1, the 
expression cos^0oK(k) vanishes as k — > 1, because it then goes as (fi;'^/2) ln(4/fi;'), where 
we have introduced the notation k' = \/\ — k}. 

The exact solution of the oscillation equation fH9|) can be written in the form 



x'y = A sin (f), (56) 

where 

^=q( , X (57) 

as follows from (l52l) and (15^ . Here, we have introduced the function Q(-2), inverse to the 
function R{(f), k), such that z = R{Q{z), k) for any real number 2; G M. This is possible be- 
cause -R(0, k) grows strictly and is continuously differentiable with respect to the variable 
on the entire real M axis when ^ k < 1, with i?(R, k) = M. When these conditions 
are satisfied, the inverse function Q{z) exists and is also a single- valued, strictly growing, 
continuously differentiable function on R, so that Q(M) = R. The derivative dQ/dz is 
never equal to zero or infinity, since this is true for the partial derivative dR{(f), K,)/d(f) 
for all C values of interest to us in the interval 1 ^ C < 00 (this corresponds to any 
physically possible value of the particle Lorentz factor 7niax)- 

Formulas (156|) and (1571) exhaust the problem of oscillations of a charged particle 
about the force-free surface in the absence of radiative energy losses. We do not require 
such an analysis of the function Q{z); let us note just one of its properties. As follows 
from the analogous properties of the elliptical integrals ( 153|) . the function R{4>, k) satisfies 
the relation 

R{7in ±4>,k) = 2nK{k) ± i?(0, k) 
for an arbitrary whole number n G Z. Thanks to this, it is true for the function Q{z) that 

Q{z + 2nR{k)) = Q{z) + Tin. 

If we pass from Q{z) to (f)(t) using fl57j) . the equality (f){t + nT/2) = (pit) + vrn follows 
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for (f){t), where T is determined by fl55|) . Hence, we can again convince ourselves that a 
charge particle will indeed undergo oscillations with the period T. For this, it is sufficient 
to take any even n and substitute the resulting relation for (j){t) into f l56p . 

Asymptotically, the functions Q{z) for k close to zero and to unity have the form 

Qiz) = \ (58) 
[ arcsin(z — 2h{z)) + nh^z), k — )■ 1, 

where the integral values of the function h{z) are determined by the formula 

z + 1 



hiz) 

Here, \_y\ denotes the integer part of the real number y. No constraints are imposed on z 
in dSS]). 

In the case of non-relativistic oscillations, we must use the asymptotic of the function 
Q{z) when k = 0. As we can see from ( 157|) . z = ujt/2 when C ~ 1. We then imme- 
diately obtain the equation x'l = Asm{ut) for harmonic oscillations with frequency u. 
In the ultra-relativistic case, C ^ 1, which corresponds to /t — 1, and substituting the 
corresponding asymptotic of Q{z) ( l58|) into (156|) yields 

x[ = {z - 2h{z)) , (59) 

where, as follows from fl57|l . z = t/A. We can see that this expression is simply another 
form of ( 150|) obtained above based on qualitative reasoning. Relations ( 150|) and (159|) are 
fully equivalent. 



5. CONCLUSION 

The dynamics of the motion of electrons and positrons in the inner vacuum magnetosphere 
of a neutron star can be represented as follows. A charged particle created far from the 
force- free surface will reach relativistic speeds over a time Trei ~ 10^^'' — 10^^^ s (|T0|) . 
and will make a transition to a quasi-stationary motion regime over a time r^j ~ 10~^ s 
f|T5l) . having traversed a distance of the order of 10 — 100 cm. The Lorentz factor is then 
7o ~ 10^ — 10*^ (|T^ . and is fully determined by the balance between the power obtained 
from the accelerating electric field, E\\/Bcr ~ 10^^ — 10^^ (ITT]) (in dimensionless form), 
and the intensity of the curvature radiation. The particle moves essentially along the 
magnetic-field line ( 127|) . since the electric drift velocity is of order Ve/c ~ 10~^ ( 12^ and 
the centrifugal drift velocity (12^ . is even smaller. The radius of curvature p and 
longitudinal electric field E\\ slowly change along the particle's trajectory, leading to an 
adjustment in the particle Lorentz factor 70. The time for this readjustment is fairly small 
(to ~ 10^^*^ — 10^^ s fl20|) ). and the particle traverses a distance of order 1 cm — 100 m 
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over this time, whose upper hmit is achieved near the force-free surface. This distance 
is appreciably less than the radius of the star, so that the Lorentz factor is essentially 
determined by coordinates of the particle. As a charged particle approaches the force- 
free surface, the quasi-stationary condition is disrupted. As a consequence, the particle 
passes through the force-free surface and begins to undergo adiabatic, ultra-relativistic 
oscillations. These oscillations decay due to radiative energy losses, while their frequency 
grows. When the amplitude Inro ~ 1 cm (j46j) is achieved, the oscillations become non- 
relativistic and harmonic, at the frequency z/ ~ 1 — 10 GHz (H5l) . Since the regular 
component of the particle's velocity lies in the plane tangent to the force-free surface, the 
particle undergoes a regular drift motion along the force-free surface, simultaneously with 
its oscillatory motion. 
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Cross sections of the force-free surface by the plane ip — ip^ = {0, tt} (left) and the plane 
^ — = {7''/2, 37r/2} (right) for angles between the magnetic and rotational axes of 0, 
7r/6, 7r/3, 7r/2. (top to bottom). The signs of the accumulating charges are indicated. 
For comparison, we also show the cross section poj — 0. 
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